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The influence of tlie liadronic tttt interaction on tlie wave functions of pionium 
(the Coulomb bound vr+vr" system) at distances < 10 fm is calculated for s-states 
with principal quantum numbers n < A. Hadronic tttt potentials are used that 
reproduce the scattering phase shifts of two-loop chiral perturbation theory. The 
pionium wave functions ipnif) are obtained by integrating the coupled Schrodinger 
equations for the (7r"'"7r~, vr^vr^) system. We find that, for r < 10 fm, n^^'^%jjn{r) is 
practically independent of n. From this we conclude that the production rates of the 
s-states of pionium are proportional to n~^, a result needed for the interpretation 
of the DIRAC experiment currently running at CERN. 
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The DIRAC experiment currently in progress at CERN[1], with the aim of 
measuring the Hfctimc ri of the Is-state of pionium with an accuracy of 10 % , 
requires for the extraction of the hfetime an accurate knowledge of the production 
probability Wn of the ns-states as a function of n. A proton beam of 24 GeV/c 
is incident on a target (for example, titanium, nickel and platinum). Among the 
reaction products are many free tt+tt" pairs and a much smaller number of tt+tt" 
pairs bound into pionium, whose production occurs almost entirely in ns-statcs. 
These pionium atoms may decay into with a lifetime r„, or they may break 

up because of their interactions with the target atoms. The effect of the decay into 
77 can be neglected here. 

The TT'^TT" pairs arising from the breakup of atoms nearly all have a relative 
momentum < 3 MeV/c in their cm. system. Observed pairs with a larger relative 
momentum are therefore almost all free pairs. This kinematical characteristic can 
be exploited to obtain experimentally for a given target the number of free pairs 
and the number of breakup pairs. By observing the number of breakup pairs from 
targets with different characteristics it is possible to deduce the lifetime ti. 

This method of measuring ri requires modelling of the processes of atom for- 
mation, decay (into n^ii^) and breakup. For this modelling it is necessary to know 
how r„ and Wn vary with n. To a high degree of accuracy, r„ = n^Ti, independent 
of the details of the hadronic interaction. This result is easily obtained using the 
formalism of Ref.[2], which describes the decay process without the use of wave 
functions. The behaviour of Wn as a function of n is the subject of the present 
paper. 

For a collision process with pions among the final state particles, let F(P, q) 
be the probability amplitude for the production of a free tt+tt^ pair with total 
momentum P and relative momentum q of the pair in its cm. frame. The formation 
of F involves an integration over the other variables that characterise the final 
state particles and a sum over all final states containing at least one tt"*" and one 
7r~. Let F(P,x) be the Fourier transform of -F(P, q) with respect to q. Since the 
production of the pair is a hadronic process, F(P,x) is appreciably different from 
zero for |x|= r < 10 fm at most. 

According to Nemenov[3], where further details can be found, the probability 
amplitude for the production of a pionium atom in an ns-state with momentum P 



where '0n(x) is the pionium wave function in the rest frame of the atom. In fact 
■0„(x) is real and depends only on |x|= r. The probability for the production of a 
pionium atom with momentum P and principal quantum number n is then 
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In the following we will show that, to a high degree of accuracy, 



i>n{r) = n ^^'^i>i{r) for r < lOfm. (2) 

In view of our earlier observation on the behaviour of F(P, x) as a function of r, it 
follows from Eqs.(l) and (2) that 

Wn = n~^wi (3) 

to an even higher accuracy. 

For a pure (i.e. point charge) Coulomb tt+tt" interaction we know that, for small 

r. 

Since the Bohr radius B for pionium is 387.5 fm, it follows that Eq.(2) holds in this 
case with an error of at most a few parts in 10^. Adding the hadronic potential 
and taking account of the finite extension of the charges results in a modification 
of ipn{i^) for r< 10 fm. Amirkhanov et al.[4] have made calculations of the modified 
ipni'f')- but they ignore the possibility of the decay of the atom into and treat 
the problem as a single- channel one. They do not modify the Coulomb potential at 
small distances and they use a simple Yukawa-type 7r"^7r~ hadronic potential. They 
conclude that Eq.(2) continues to hold to a high level of accuracy. 

Our work improves on Ref. [4] in three important respects. We have included 
the Coulomb potential for extended charges, exactly as in Ref. [5] (last paragraph 
of Section 4). We have used more realistic hadronic potentials and we have taken 
account of the presence of the tt^tt*^ channel and used coupled Schrodinger equations 
to model the (tt+tt", 7r°7r°) system. For the s-wave and for 2/io <W< 2//c {W is 
the total energy of the tttt system in its cm. frame and and are the masses 
of the charged and neutral pions) these may be written in the matrix form 

(l2^ + Q^-MfV(r))u(r) = 0, (4) 

where 

V(r) = V"'"(r) + V^{r) . (6) 

Thus qo is the cm. momentum in the 7r°7r° channel. The electromagnetic potential 
matrix is 

/ yem q \ 
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Here 



where V^'^ is the Coulomb potential between extended charges; there is no need 
to take account of the very small relativistic modification of the Coulomb potential 
or of vacuum polarisation. The hadronic potential matrix is assumed to be isospin 
invariant, in accordance with the arguments of Gasser and Leutwyler [6,7]. Thus 

where is the hadronic potential for / = and total isospin /. 
The diagonal matrix f has the nonzero components 

jc — 77} — , JO — 77} — ■ (,y; 

Its inclusion in the Schrodinger equation using a potential has been extensively 
discussed in Refs.[8] and [9]. Since we will work very close to the 7r~^7r~ threshold, 
where the pionium states are located, it is sufficient to take 

The column vector u(r) is 

where 

Uc{r) = V4nrip{r) , Uo{r) = rRo{r) , (11) 

ilj{r) being the tt+tt^ wave function and Ro{r) the tt^ti^ radial wave function. For 
each value of < 2/ic, there is a unique solution of Eq.(4) which satisfies 

u(0)=0, (12) 

as well as the asymptotic behaviour for r — > oo 

ryrri 

u,(r) ^ A,r^e-^^\ C^^, (13) 

Uo{r) = Ao sm{qor + 6oo) , Aq > , (14) 
and in addition the normalisation condition 
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dr{uc{r)f ^ 1 . (15) 



Here Sqo is the phase shift for elastic n^n^ scattering and depends on 1^. If 
W goes through energies Wn that are close to the energies of the pure Coulomb 



bound states, Sqq goes through an odd multiple of 90°. These energies are the bound 
state pionium energies with the influence of the hadronic tttt interaction taken into 
account; the functions V-'n('") at the energies Wn are the modified wave functions in 
which we are interested and which we have calculated for n < 4. 

This completes the details of the calculation of the wave functions i/Jnir), except 
for the hadronic potentials V^{r). We have constructed potentials that are inde- 
pendent of both the energy W and the hadronic mass ^, which reproduce very well 
the phase shifts given by two-loop chiral perturbation theory [10,11] for /i — /iq and 
/J, = He via the Schrodinger equations 

+ ql - l^c fc V^{r)Wir, = , (16) 

+ Qo - /^o /o V'ir))u'{r, /^o) = , (17) 
with 1 — 0,2. The factors fc and /o are given in Eq.(9) and by 

ql^ — -^J,l , W>2iic. 

The potentials V^{r) are of the same analytical form as those used for the 
calculation of the electromagnetic corrections in low energy pion-nucleon scattering 
[8,9]. In addition to a range parameter that was fixed at 1.5 fm, there are three 
further parameters that were varied to obtain the best fit. Full details of this 
parametric potential are given in Ref. [12]. The phase shifts were generated by 
integrating the regular solutions of Eqs.(16) and (17) outwards from the origin. The 
potentials F°(r) and V'^{r) are plotted in Fig.l. We thus have captured part of the 
dynamical behaviour of chiral perturbation theory for tttt scattering by means of a 
simple model that uses potentials independent of the energy and the hadronic mass. 
This makes possible a reliable calculation of the pionium ns-state wave functions 
and therefore of the production rates Wn as a function of n. 

To show that n^^'^'ipnir) ^ i^i{r) for r < 10 fm for a range of values of n, we 
have calculated the quantity 

dnir)^n^/^iPn{r)/A{r)-l, r<10fm, 2,3,4. (18) 

Fig. 2 shows the result of this calculation. The values at r = arc ^2(0) = —0.0014, 
dsiO) = —0.0021, 0^2(0) = —0.0022. As r increases, dn{r) changes quickly at first 
and then becomes almost constant. The value of d„(r) increases in magnitude as 
n increases, but Fig. 2 shows that Eq.(2) holds for n = 2, 3, 4 to better than 3 parts 
in 10^. The corresponding results for n > 4 are not important, since only 2 % 
of all pionium atoms are produced in states with these higher principal quantum 
numbers. 

The results of Ref. [4] have the same sign as ours and show the same rapid change 
near r = 0. At r = 8 f m , Ref. [4] gives ^2 = -0.0012, = -0.0010 compared 



with our values of -0.0020 and -0.0026 respectively. Our more refined calculation 
yields results that are larger in magnitude but Eq.(3) still holds with an accuracy 
completely sufficient for use in modelling the DIRAC experiment [f]. 
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Figure 1: The potentials V^{r) and V^{r) that give the best fit to the two- loop 
CHPT phase shifts of Ref.[10,ll] for hadronic masses pLc and /Xq- 
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Figure 2: The functions rf„(r) defined in Eq.(18) . 



